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Intro PhysX' nVIDIA 


• 1999 Ph. D. ETH Zurich: Polymer simulation 

# 1999-2001 Post doc MIT: Simulation in CG 

% 2002 Co-founder NovodeX (physics middleware) 
% 2004-2008 Head of research Ageia (SDK features) 

% 2008 Research lead PhysX SDK at Nvidia 

m Zurich office: R&D PhysX SDK 
<i> currently Height Field Fluids 
<§ New Ideas 
# Lessons learned so far 
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Hello World! 


PhysX nviDiA 


% Use two arrays float u [N,M] , v[N,M] 

% Initialize u[i,j] with interesting function 
% Initialze v[i, j]=0 

demo 


loop 




+=(u[i-l, j 

+ u [i+1 , j ] + u[i,j-l] + u[i,j+l])/4 - u[i,j] 

v[i,j] 

*= 0.99 


u[i, j] 

+= V[i,j] 


visualize (u [ ] ) 
endloop 



Clamp at boundary e.g. def. u[-l,j] = u[0,j] 
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Motivation Physx nviDiA 


That’s it: Thank you for your attention! 

<§> The magic (physics) behind it 
# Object interaction 
% Boundary conditions 
% More complex domains 
% Barking waves 

(§> Implement it yourself (use CUDA!) 

% Soon: 

<§ Use the PhysX SDK: scene->createHeightf ieldWater ( . . ) 

• Drop HF water into your game 
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Outline 

% Introduction 
S Fluids are cool! 

# From offline to real-time 

% Physics is not hard 

• f = ma 
II Simulation 

% Height field water 

% Simulation of columns 
<§ Boundaries 
S Object interaction 
% Horizontal motion 
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<$u 

PhysX nviDiA 



Off-line Fluid Simulation 


PhysX nviDiA 


% State of the art is impressive! 

% Google “Robert Bridson”, “Ron Fedkiw”, 
“James O’Brien”, ... 


Liquids 
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Example 


PhysX nviDiA 



Guendelman et al. (SIGGRAPH 2005) 
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Off-line Simulation Times 



# Typical setup 

# Grid size 51 2 3 cells 

# Linear system with 
>100 million unknowns! 

# Level sets on even finer grids 

# Raytracing 

(reflection / refraction / caustics) 



# Photorealistic results 


#10 seconds - 50 minutes per frame! 


# Recently real-time on a GPU (Crane et al.) 

# Coarser grid, simplified physics and rendering 
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Game Requirements 

PhysX nviDiA 

% CHEAP TO COMPUTE! 


# Small fraction of the 15 ms of a frame 


49 Stable even in non-physical settings 


41 Kinematic, fast moving objects / characters 


49 Low memory consumption 


% Challenge: 


41 Get as close as possible to off-line results 
9 Meet all these constraints! 
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Solutions 

4§> Resolution reduction 
4§ Blobby and coarse look 
9 Details disappear 

49 Use real-time tricks! 

9 Reduction dimension 3d -> 2d 
9 Physics low-res, appearance hi-res (shader effects) 
41 Simulate active and visible regions only (sleeping) 
9 Level of detail (LOD) 


PhysX nviDiA 
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Approaches 


PhysX nviDiA 


• Procedural Water 

# Unbounded surfaces, oceans 


<§ Particle Systems 

# Splashing, spray, puddles, smoke 


H Height field Fluids 

# Ponds, lakes, rivers 



Procedural Animation 



% Simulate the effect, not the cause 

% [Fournier86], [Hinsinger02], [Bridson07], [Yuksel07] 


# Pros 

% High level of control 

# No limits to creativity 

(e.g. superimpose sine waves) 

# Cons 

% Modeling fluid-scene interaction 
more difficult 

# Less realistic 



itfQGbt demo 
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Particle Based Fluids 


PhysX nviDiA 


•' Fluid is represented as a set of particles 

# [Monaghan92], [Premoze03], [Muller03] 

©Pros 

• Particle simulation is fast and quite simple 

# Spray, splashing, small puddles, bl 
runnels, debris 

# Cons 

% Surface tracking difficult 

• Screen space meshes [Muller07] 

% Not suited for lakes, rivers, oceans 
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Physics for simulations 
is not that hard 
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Useful Equations (PDEs) 


PhysX nviDiA 


% Fluid (gas / liquid): 

p(v, + v-Vv)— -Vp+f 


SH V V ' V V ' 

m a = f 

% Elastic solid: 

p U„= V-cr 5 (u)+f 


m a = f 



Newton’s Second Law 




change of velocity 
per unit time 


force divided by mass 


Simulation Loop (explicit Euler integration) 

# compute force 

(based on actual positions and velocities) 

% velocity += force/mass time step 
% position += velocity time step 
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Example: Particle 

PhysX nviDiA 

9 Need to store position and velocity! 

v, -At 

V 1 + 


m 

\v~\t ( 

|L. V 2 




f a 


Vj'AN^ 
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Other Particle Forces 


PhysX nviDiA 


Smoke, debris 



Liquid (e.g. SPH) 


9 — C 


particle ^ game level! 


Neighbor search! [Teschner03] 
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PhysX nviDiA 


Height Field Fluids 
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Height Field Water 


PhysX nviDiA 


continuous 


discrete 


<§> Represent fluid surface as a 

Hi 2d function u(x,y) 

# 2d array u[i,j] 

(§> Pro: Reduction from 3d to 2d 

# Cons: One value per (x,y) -» no breaking waves 






Position and Velocity 


PhysX nviDiA 


•; Position = column height 

$ Velocity = change of column height per time unit 
% (h = column width) 

\v = u t 

u I 
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Force 



# Force causes waves to travel at speed c 
in all directions 



demo 
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Constant Slope 


PhysX nviDiA 



# Constant slope no velocity change -» no force 
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pos curvature up acceleration neg curvature down acceleration 


->• force proportional to curvature! 
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Wave Equation 


PhysX nviDiA 


# Newton’s Law (continuous): 

«« =f/m = ku xx /m 

# Assume m constant : Can replace k/m by c 2 

u tt = c2u xx Id wave equation [Jeffrey02] 

# Solution: 

u(x,t) =f(x + ct)+g(x- ci) for any fg 

# Constant c is the speed at which waves travel 
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u x [i-%] = (u [i] -u [ i— 1 ] ) /h u xx [i] = (u lt [i+y 2 ] -u x [i-y 2 ] ) /h 

u x [ i+%] = (u [i+1 ] -u [ i] ) /h = (u [ i+1 ] -2u [ i ] +u [ i- 1 ] ) /h 2 


Sid curvature: 

(u [i+1] +u [i-1] -2u [i] ) /h 2 

S 2d curvature: 

(u [i+1 ,j]+u[i-l,j] +u [i , j+1] +u [i , j-1] -4u[i, j] ) /h 2 
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Column Simulation Step 


PhysX nviDiA 


• Id simulation 



% 2d simulation 



Some Remarks 

PhysX nviDiA 

% Information can only propagate 1 cell per time step 

% Upper limit for choice of wave speed 


% c<h /At 


% or upper limit for choice of time step 


% At < h/c (Courant-Friedrichs-Lewy (CFL) condition) 

% Boundary conditions needed 


% Clamp yields wave reflection 
9; Wrap yields periodic propagation 
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Numerical Explosions 


PhysX nviDiA 


•; Explicit Euler step x += v-At 


% Assumes velocity constant during a time step 
% Particle sample: % Column sample 



GameDevelopers 

Conference 

08 


Damping vs. Clamping 



• Damping force (physically “correct”) 

• f = -tv 

• But how to choose k? No stability guarantees 


Scaling (unphysical, more direct control) 

v = sv // s < 1 , smaller time step -> stronger effect 

® Clamping (unphysical, direct control) 

offset = (u[i+l, j]+u[i-l , j]+u[i, j+1] +u [i, j-1] ) /4 - u[i,j] 
maxOffset = maxSlope * h // independence of resolution h: 

if (offset > maxOffset) u[i,j] += offset - maxOffset 
if (offset < -maxOffset) u[i,j] += offset + maxOffset 
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Object Interaction 


PhysX nviDiA 


# Object water 

# Object pushes columns beneath it down 
% Add the removed water in the vicinity! 



% Water ->• object 

% Archimedes' principle 

# Each column below the object applies force 
f = -Au h 2 p g to body at its location 

% Au is the height replaced by the body, 
p water density, g gravity 
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Fully Immersed Bodies 


PhysX nviDiA. 


% Body below water surface 
% Hole appears above the body 
# Non-physical 


• See story of divided sea 
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Solution 


PhysX nviDiA 


•; New state variable 

d Each column stores the part 

of u [ i , j ] currently replaced by solids 


% At each time step: 

« u[i,j] is not modified directly 

m 

is distributed as water u 
to the neighboring columns 
d In case of a negative difference 
water is removed 



GameDevelopers 

Conference 

08 



GameDevelopei 

Conference 





Water Rendering 


PhysX nviDiA 


% Reflection 


T 



cube map 


% Refraction 


n 


<- 


Scene 

below surface 
rendered to 
texture 


% Caustics: Cheating - Animated texture 
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Unbounded domain 


PhysX nviDiA 


% Simulations in tank not practical for games 
% Solution: Hybrid Simulation 



3d Euler 



Height field 


Height field 

Procedural 


[Thuerey06] 

% No reflections or wrap on boundary! 
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3d-2d Simulation 

PhysX nviDiA 

9 N. Thuerey et al., Animation of Open Water Phenomena with 
coupled Shallow Water and Free Surface Simulations, SCA 06 



r* 

Wyi 

3D domain highlighted 
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Open Boundaries PhysX' nvioiA 


% Let waves with speed c pass freely 
% Use ghost column on border with height memory g 



C 


% Temporal change of g proportional to border slope 

(g new - g)/At = -C (g new - u)/h 

• Update rule: 

g neM = (c-At-u + gh) / (h + c-At) 
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Demos 


PhysX nviDiA 







Horizontal Motion 


PhysX nviDiA 


%■ We only work with vertical velocities 
% Horizontal velocity for dragging floating objects 


v[i] 



v h [i] = (v [i] -v [i-1] ) h + (v [i+1] -v [i] ) h 


% Better model: Shallow Water Equations 
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Shallow Water Equations 

PhysX nviDiA. 

DulDt = -uVv 
D\/Dt = -g V(u+b) + a ext 


• Oversimplified: 

9 Compressible Eulerian fluid simulation 
on a 2d staggered horizontal grid 
C Density interpreted as height 


% No further details at this time © 
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Water Loss 


PhysX nviDiA 


% Big problem in Euler / level set simulations 
% Less problematic in height field fluids 

% Simple, bullet proof solution: 

# V = E i(j u [i , j ] 

® V new = 2 i,j U new [i, j] 

# Distribute v-v new evenly 

# Within connected regions! 
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•; Height fields cannot capture breaking waves 
% Identify steep wave fronts 
% Construct and track line along front 
• Emit patch of connected particles 

% Generate mesh with given thickness 
for rendering (plus particles for foam) 
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Breaking Waves 


PhysX nviDiA 






Surfing... 

Simulation Resolution: 
Avg. Simulation FPS: 


200x100 

40.6 
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PhysX' nviDiA 


Questions? 


Slides available soon at www.MatthiasMueller.info 
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